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. Abstract 

O 

We present the coupled nonlinear integral equations (NLIE) governing the finite 
size effects in N=l super sine-Gordon model for the vacuum as well as for the excited 
states. Their infrared limit correctly yields the scattering data of the super sine- 
Gordon S-matrix conjectured by Ahn Ultraviolet analysis is in agreement with 
the expected conformal data of c = 3/2 CFT. Conformal perturbation theory further 
corroborates this result. 



1 Introduction 

Finite Size Effects in two dimensional integrable QFT can be studied exactly by means of 
sets of coupled nonlinear integral equations that can be thought as continuum versions of 
Bethe Ansatz equations. In particular, the equation introduced for the Sine-Gordon (SG) 
model by Destri and De Vega [HE] through a light-cone lattice regularization and basically 
equivalent to a similar one found by Kliimper, Batchelor and Pearce [SI El for the spin 1/2 
XXZ chain, once extended to excited states E], has proven to be a very important 
tool in establishing and testing the correct connection between the Factorized Scattering 
description and the Perturbed CFT formulation of SG model IH1 EH- Restricted SG 
models can be put in correspondence with minimal CFT perturbed by their 0i3 operators 
in a similar way [Tfl] . 

It would be a result of principal importance to be able to formulate such non-linear 
integral equations (NLIE) for many other integrable QFTs. In this paper we report the 
construction of NLIEs for the N=l Super Sine-Gordon (SSG) model. The importance of 
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this generalization lies in the manifest supersymmetry of this integrable model and the 
link it has with applications in various areas of physics ranging from condensed matter to 
strings. 

To get the NLIEs for SSG model, we proceed by analogy with the pure SG case. There, 
the NLIE was obtained through a light-cone lattice construction of a 6-vertex model, 
which is known to be equivalent, from the Bethe Ansatz point of view, to an alternating 
inhomogeneous spin 1/2 XXZ chain. In the continuum limit this lattice model is shown 
to give the equations of motion of SG (or equivalently massive Thirring) model [TTlfH] . 
It is known since a long time that the 6-vertex model on a square lattice, or equivalently 
spin 1/2 homogeneous XXZ chain, renormalizes on the continuum to a free massless boson, 
i.e. c = 1 CFT [THj. The analog model with spin 1, the integrable spin 1 homogeneous 
integrable XXZ chain is equivalent to the square lattice 19- vertex model and renormalizes 
to a c = 3/2 boson + Majorana fermion theory which is explicitly N=l supersymmetric. 

It is therefore natural to conjecture that the full SSG model can be constructed starting 
from a light-cone lattice 19-vertex model, or equivalently from a spin 1 alternating inho- 
mogeneous integrable XXZ chain [2B]. A direct deduction of the equations of motion, like 
in [TT] or [12] is too cumbersome and beyond the scope of the present paper. 

The Bethe Ansatz for spin 1 integrable XXZ chain has been written long time ago [16j 
and can be generalized with the introduction of an alternating inhomogeneity trivially. 
In the homogeneous case, a coupled set of two NLIE's equivalent to it on the lattice has 
been deduced, by use of transfer matrix recurrence relations and T-Q formalism in ref . [Tl] . 
including excited states. Their generalization to the inhomogeneous case is also trivial. On 
the other hand a system of two NLIE was conjectured on the continuum by C. Dunning (2^1 
for the vacuum state and checked to give the expected ultraviolet central charge c = 3/2. 
Further analysis of Finite Size Effects in SSG model was carried out in [18J. 

In section 01 we recall the setup of the model in Lagrangian formulation, as a perturbed 
CFT and a factorized scattering theory. In section 03 we present the alternating inhomoge- 
neous integrable spin 1 XXZ chain and its Bethe ansatz. The NLIE is presented in section 
El for the spin chain on the lattice and, by performing the scaling limit, in the continuum. 
Section El is devoted to the infrared (IR) analysis reconstructing the supersymmetric soli- 
tonic S-matrix. From this section on we limit ourselves mostly to the repulsive regime. The 
interesting features of the attractive regime, as well as the implications for superminimal 
model perturbed by $i 3 are left for a future investigation. In section|Hlthe ultraviolet (UV) 
limit of the NLIE is computed and shown to reproduce data compatible with the expected 
c = 3/2 CFT. A numerical comparison with predictions of conformal perturbation theory 
is presented in sectional Finally we draw our conclusions and perspectives for future work 
in section [HI 
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2 The super sine-Gordon (SSG) model 

The SSG model can be defined on a cylinder of circumference L (thought as spatial direc- 
tion, with infinite time direction) by the action 

Assg = f dt f o dx J Md ® ( D ®D<5> + n cos p$) 
where (3 is the coupling, fi a bare mass and the superfield $ expands in components as 

t, 9, 6) = tp(x, t) + $i/j(x, t) + {hj){x, t) + MF(x, t) 

The action in components, once the auxiliary field F is eliminated by the Euler-Lagrange 
equations of motion, reads 

Assg = f dt f o dx {k^^f + i$'td v ti> + V^>i> cos ^ + ^ cos fiipj (1) 

This model is known to be integrable since a long time and evidently possesses N — 1 
supersymmetry. Its ultraviolet (UV) limit is given by a theory of a free massless boson 
Q(x, t) = (f>(z) + <j)(z) = (here z = x + it, z — x — it) compactified on a circle of radius 

R = h£ (i.e. with quasi periodic boundary conditions Q(x + L,t) = Q(x,t) + 2irmR 
at fixed time, m being called the winding number) and a free massless Majorana fermion 
with antiperiodic (NS = Neveu-Schwarz) or periodic (R = Ramond) boundary conditions 
^(x + L, t) = ±tp(x,t). This CFT obviously has central charge c = | and a (U(l) x Z 2 )l x 
(£7(1) x 1^2)r symmetry and shows N = 1 supersymmetry. 

In the NS sector the primary fields are given by vertex operators 

Vi%\z,z)=^ r : e *[(ft+ 1 * l )*M+(ft- H iW«>] : 

where r, f G {0, 1}, ipo = 1, ip\ = ip of conformal dimensions 

a(^-) = I (Hl + !^V + !1 a(^) = -(-- — V + f 

n ' m 2\i? + 2 j + 2 ' n ' m 2 Vi? 2^2 

In the R sector the primary fields are realized by bosonic vertex operators multiplied 
by the Ising spin field 

R n , m (z,z) = a(z,z) : e'Ki+^Kfc-W*)] : 
of conformal dimensions 

1 fn mR\ 2 1 T 1 fn mR\ 2 1 



2\R 2 J 16 ' 2 \i? 2 / 16 

All the other secondary states in the theory can be obtained by applying the Heisenber^ 
algebra and the fermion algebra modes to the primaries. See [TH] for details. 
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The main object encoding the UV operator content of the theory is the modular in- 
variant partition function of the c = | CFT, obtained in |T2j. It reads 

Z ( R ) = I (XoXi/2 + X1/2X0) 

/ 1 ,9 1 | 2 \ 1 ,2 A (0 ' 0) A (0 ' 0) 

+ (IXoh + IX1/2I 2 ) 2^ +IX1/I6I 2^ [ g n >™ g n > m 

nSZ,mG2Z 2n-mS2Z+lJ 

where 77(g) is the Dedekind function and g = e 27nT , r being the modular parameter. We see 
from the partition function that in the NS sector the quantum numbers n, m, r, f describing 
the operator content can take the values 

n G Z + - m G 2Z + 1 <^ n _ 1 
2 [ r = r = 1 

or 

[ r = r = 

n g Z m g 2Z 



r = r = 1 

while in the R sector they are constrained to take the values 

n G Z m G 2Z + 1 

or 

n G Z + - m G 2Z 

The SSG model can be seen as a perturbation of this c = 3/2 CFT by a relevant 
operator 

% ert = V2n {v}^ + Vilj) = v^vr : ^ cos ^ : (2) 
of conformal dimension A pert = A pert — | + The other perturbing field 

necessary in the classical action to ensure supersymmetry, can be shown to give a contribu- 
tion to quantum correlators that goes to zero with the cutoff [21] and can be neglected in 
a quantum formulation. See |18| for a more detailed discussion of these two perturbation 
terms. 

The operator (j2j) is relevant (A pert < 1) if the SSG coupling is in the range < (3 2 < 16n. 
At exactly f3 2 = 16ir it becomes marginally relevant and describes a perturbation of the 
WZW S77(2) fe=2 model by its J a J a operator. 

Opposite to this UV description of SSG as a perturbed CFT, there exists a factorized 
scattering theory that has been conjectured by C. Ahn [TH] to describe the asymptotic 
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particle properties of SSG model. The S-matrix is written in detail in the original paper 
[T5] and the vacuum structure is carefully analyzed in [THJ. While not repeating here those 
data, we just recall that the S-matrix has a tensor product form 

S(e) = SRsos(2)(0)®S SG (9,l3sG) 

where the RSOS part is the S-matrix of the tricritical Ising model perturbed by ^13 and 
describes the supersymmetric structure, while the SG part is the usual solitonic S-matrix of 
pure sine-Gordon model, at 0g G = , describing the bosonic degrees of freedom. The 

vacuum structure shows, similarly to SG, infinitely many degenerate vacua periodically 
repeating. The additional RSOS structure amounts to the fact that in each period there 
are three degenerate vacua, labeled by {0, |, 1} and connected by the adjacency of the 
A3 Dynkin diagram. The elementary excitations are supersymmetric solitons (kinks) that 
interpolate these vacua [IHj. They are the only asymptotic particles of the model if ^ < 
P 2 < 167r (repulsive regime). For [3 2 < bound states (supersymmetric breathers) 
appear (attractive regime). 

3 The inhomogeneous spin 1 XXZ chain and its Bethe 
ansatz 

As said in the introduction, we assume, by analogy with the SG case, that the SSG model 
can be constructed by taking the appropriate continuum limit of an integrable lattice model: 
the light-cone lattice 19- vertex model. For a description of this lattice model including its 
Boltzmann weights, see [17] . Actually, there are infinitely many possible lattice models 
providing SSG as a continuum limit. Among them, it is convenient for our purposes to 
choose one that preserves integrability on the lattice, so that the Bethe ansatz machinery 
can be used to diagonalize its transfer matrix and get the energy levels and in principle 
other measurable quantities. 

The square lattice 19-vertex model, in the limit of continuous time direction while 
keeping the space direction discretized , is related to the integrable homogeneous spin 1 
XXZ chain [21] with anisotropy 7. Its continuum limit is a system of a free massless boson 
compactified on a radius R = ^ti/(tt — 27). The light-cone lattice version of the same 
model is instead equivalent to a spin 1 XXZ chain, with N sites labeled by n — 1, ...,N, 
lattice spacing a, anisotropy 7 and inhomogeneity A n = (— l) n+1 6. The relation between 
the compactification radius R and the SSG coupling {3 implies that the anisotropy 7 is 
related to j3 through ^ = A pert ((3). It is often conveniently parametrized by 7 = 
Throughout this paper, we will mainly focus on the range 7 < |, i.e. p > 1. This 
corresponds to the repulsive regime of SSG model. 

We start from the description of the homogeneous system. At each site n of the chain 
a 3-dimensional space V n ~ C 3 is defined, on which the spin operators of the Hamiltonian 
act. Let Rn,m(0) be the R matrix acting non trivially on 14® V m . Its explicit form can be 
found in [21]. 
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The homogeneous transfer matrix T 2 {9) is defined by the trace of the product of R over 
the auxiliary space V a ~ C 3 

T 2 {9) :=Tv Va R a , 1 (6)---R a , N (6). 

Periodic boundary conditions are imposed here. 

The hamiltonian of the associated spin 1 XXZ chain is identified with H = logT 2 (0). 

TV 

n =H K ~ + cos2 7 « - «) 2 ) - (2cos 7 - 1)(^< + - 4sin 2 7 (^) 2 ) 

n=l 

(3) 

where we write a n = S n ■ S n+ i = + cr*. 

Next we consider the inhomogeneous generalization (TJ [2j [2H| , (N £ 2Z) 

T2(0\{Xj}) := Tw^a^e - X,) ■ ■■R a , N (9 - \ N ). 

By following the light-cone approach, we adopt a particular alternating choice of the in- 
homogeneities A n = (— l) n 6 and write T 2 (9,Q) = T 2 {9\{{— l) n 6}). The local hamiltonian 
containing both the two-body and the three-body interaction is then defined by 

H=--^\ogT 2 (9)\ e=e 
7 at) 

The inhomogeneity introduces the interaction between the left-going frame and the right- 
going one. 

Conveniently, we introduce an auxiliary transfer matrix which is defined in a similar 
manner to T 2 , replacing the spin 1 x spin 1 R by the spin 1 x spin | R matrix. A twist u 
at the boundary can also easily be incorporated. Periodic boundary conditions correspond 
to uj = 0, antiperiodic ones to u> = ir/2. The explicit eigenvalues of Ti{9) and T 2 {9) then 
read 

me) = ^p{9-i^)9^l + ^p { e + i^ e - i ^ 



Q(9) rK ' Q(9) 



7T 



37r.Q(0 + zf) ,„ .tt, ,„ ,7r,Q(9 + if)Q(9-0 



T 2 (9,e) = e^p(9 - i-)p{9 - i— ) ^ /' + p(9 - i-)p{9 + %-) 



2 > Q{9-%\) ,y 2 ,ry 2' Q(9 + q)Q(9-i 



2 • 



2" y 2 ' Q(9 + i w 



2 



where 



/2 2 ( ^_ e)sinh iv/2 2 ( 

7T 7T 

The important function Q is given by the Bethe ansatz roots 9j,(j = 1, • • • , M); Q(z) 
tM 



p(z) = sinh^ 2 ±{z-&) sinh^ 2 ±(z + 0) 

;ior 

UT^h^z-9,). 
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Bethe Ansatz for the homogeneous case was written in (THJ . For the alternating inho- 
mogeneous case it can be generalized easily, resulting in 

/ sinh i(0j + 9 + ztt) sinh ^ - 9 + in) \ N/2 _ 2iuj smh^ - 9 k + in) 
V sinh I (6 3 + Q-in) sinh £ (0. - 9 - m) J 11 sinh 2 - fe - 

The number of roots M is related to the (conserved) total 3-rd component of the spin 
S z of the chain and the number of sites iV by 

S z = N — M 

The Bethe equations have a periodicity for 9-^9 + ^ = 9 + in(p + 2). It is convenient 
to consider as fundamental strip — ^£±21 < \ m Q < li^+Hi 

The vacuum state is realized by a maximal set of (quasi)-two-strings (see Fig. 1) 

9 j = Re9 j ± i (| + e,) 

where the deviations 6j > from the pure two-string behavior are assumed to be small 
and limited by a positive parameter e: ej < e. The distribution of these roots is denser in 
the regions near ±9, where also the 6j are minimal. We exemplify this in some plots, at 
p = \ (attractive) and p = | (repulsive) in appendix C. 

The excited states are usually characterized by roots located in positions different than 
those of the vacuum, as well as some holes. The following classification of the Bethe ansatz 
roots, other than 2-strings, is possible: 

1. inner roots : |Im^| < |, j = 1, • • • , M I 

2. close roots: | < |Im0j| < ^, j = 1, • • • , M c 

3. wide roots: ^ < |Im^| < f^, j = 1, • • • , M w 

2 

4. self conjugate roots : |Im0j| = |-, j — 1, • • • , M p . 

Besides these roots, the real zeroes of the transfer matrices T\ and T 2 also play an impor- 
tant role in the description of the excited states. The zeros of transfer matrix T 2 (9) are 
called holes because they are in 1:1 correspondence with actual holes in the quasi 2-string 
distribution. 

Once the Bethe roots are calculated as solutions to the Bethe equations, the formula 

p q(E±P) = +iw TT sinh ^7T ± 3 - - 9) 
1 = 1 sinh 1(*7T T ^ - 9) 

yields the energy level and total momentum of the corresponding state. Higher integrals 
of motion can also be expressed in a similar way. 

We remark that the generalization of the above formulation to general values of spin is 
immediate. 
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4 The NLIE for the inhomogeneous spin 1 XXZ chain 

NLIE on the lattice for a homogeneous spin 1 XXZ chain with periodic boundary conditions 
has been written in [H], to which we refer for notations and presentation of the deduction. 
The generalization to the alternating inhomogeneous lattice model is quite straightforward 
and leads, along similar lines, to the following coupled NLIE 

In b{9) = C b + iD b {9 + ie) + ig x {9) + ig b (9) + (G * In B) (9) - (G [2e] * In B) (9) 

+ (K [ -^ +e] *\nY)(9) (4) 
In y{9) = C y + g y (9) + {K^~ e] * In B){9) + (K^ +e UlnB)(9) 
where B{9) = 1 + b{9) , Y{9) = 1 + y{9) and for any function / and g 

(f*g)(x)=[dyf(x-y)g(y) , f™(x) = f{x ± irf) (5) 



-oo 



The kernel functions G and K read 



2tt 2sinh^coshf v ' 27rcosh(#)' 

— oo 

We also introduce odd primitives of kernel functions, 

e e 
X{9) = 2% J dx G(x) xk{9) = 2% J dx K{x) 



that are important in writing the source terms containing information on the excitations: 

N H N§ M S 

9b(0) = ^x(0-h J )-2'£x(0-hf)-'£ l (x(0-8 j ) + x(0-8 j )) 
j=i j=i j=i 

M c M w M sc 

j=i j=i j=i 

9i(9) = J2xK(0-hf ] ), 
g y (9) = lim g v [ 9 + i- - %t\ 

r7 — >0+ V Z / 

E " ^) - 2 E - *?) - E - »i) + Xk{9 - 8j)) 



v- 

N„ Nfj m s 



j=l j=l j=l 

M c 

^XK(9-Cj) 

3=1 



where the second determination function, for p > 1 



Xii(d) 



In sinh ( - ) — In sinh 

P. 



9 + Sign(Im#)27r 



P 



is defined by x( x ) + xi x ~ wsgn (x)). 

The driving term bulk contribution exists only in the equation for In b{9) and it reads 



D b {9) = N arctan 



sinh 9 
cosh 



The constants C& and C y are determined by the asymptotic behavior of both sides of 
the NLIE 



C b = iC b C b = n5 b + a a = lu ( 1 + - ) + Xoo (M - N+). (6) 



7T / 1 

5 6 G {0, 1} and Xoo = = 77 1 ~ ~ 

2 V p 



(7) 



N. 



± 



5, 



u; 
=F - 



,P+2 7T 

where |_^J stands for the integer part of x and 

C y = m8 y + m(S + M w + M, 



S z 



=F - 



p + 2 77 



5 y = iV_ mod 2 



(8) 



(9) 



We need one more equation for the determination of wide and self-conjugated roots. 
For the cases n < lm9 < |(p + 1): 



lna(0) = C a + z^ a (0) + (G'7 J * lnfl)(0) - {Gfj * ln5)(0) 



(10) 



where 



9a{9) 



N H Nfi M s 

E Xii(0 - hj) - 2 XiiW - hfj - E (Xii(0 ~ *i) + Xn(9 - s,)) 
j=i j=i j=i 

M c M w M sc 

E - 9) - E - " E - „ , (11) 

j=i j=i j=i 



C- a = i{2uj[ 1 + +-(N + -N.) 
P ) P 



(12) 



and the "second-second" determination is defined as in [Tfl] . 

The parameters hj, Sj, Cj, ... in the source terms, i.e. the positions of holes, close, wide 
roots, etc... can be determined recalling that all these objects annulate the functions 
B, 1 + a, Y, etc... This leads to the following quantization conditions 
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• For holes: 

- Inbihj - ie) = 2tt I h , j = 1, N H . (13) 

• A bit formally for special roots: 

- ]nb( Sj -ie) = 2tt/ Sj j = 1,...,M s . (14) 

% 

• For close roots (only for the upper part of the close pair): 

Tln6(ct-ie) = 27rJ c T j = 1, M c /2. (15) 

• For wide roots: 

t lna(wj) = 2nl w , j = 1, M w /2. (16) 

• For self-conjugated roots: 

- lnaK^')) = 27r^TCi) J = 1, M sc . (17) 

So far we have determined only the upper part of the complex pairs, but the other parts 
can be determined by simple complex conjugation. 

• For zeroes of T\ {9) : 

ylnj/ 1 (/if ) -z|)=27r/ A (i) .7 = 1,...,^. (18) 

All the above quantum numbers J^.'s are half integers. A state is then identified by a 
choice of the quantum numbers (1^. , I c ., ...). The NLIE itself can impose constraints on 
the allowed values some of these quantum numbers. 

A comment must be done to explain the emergence of special objects (holes and/or 
roots) hj and Sj,Sj in the equations above. Normally the function Z{6) = Imlnb(#) is a 
monotonically increasing function of 9. It may happen, however, that this monotonicity 
is locally violated. If it happens that a root or hole 9j has Z(9j) < 0, then we say that 
it is special. In this case the term InB = ln(l + b) inside the convolution terms goes 
off-branch, so that InB = log£> + 2ni (where we use the notation log to indicate the 
fundamental branch of the logarithm). The 2ir jump is responsible of the emergence of the 
new terms x(9 ~ s ) + x(6 ~~ s) (or 2\{9 — h s ) for holes) in the source part of the equation. 
Smoothly varying parameters in the NLIE, if an object that was not a special becomes 
now special, it is compelled to be accompanied by two new emerging holes. So the number 
N e ff = N H — 2N S is conserved by smoothly varying parameters. 
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A careful analysis of the 6 — > ±00 asymptotics of the NLIE shows that the numbers 
Nh, Mc, Mw, ■■■ etc... are not arbitrary, but connected to the total 3-rd component of spin 
S z through a set of Diophantine equations, the so called counting equations 

N H - 2N S = M c + 2(S Z + M w + M sc ) - N+ - AL, (19) 



N l -2N* i = S z + M w + M s 



where M c is the number of those "close roots" whose imaginary part are either in {ti/2, tt} 
or {— n, — 7r/2}. There is an additional parity constraint on N\ 



M, 



R 



p+2 



to 

71 



s z 



p 



IT 



■ (20) 



S, 



P 



7T 



s z 



P + 2 



7T 



even. 



(21) 



We are interested in the scaling limit N — > 00, a — > 0, while the circumference of the 
cylinder L = Na remains finite. This limit defines the theory on a continuum cylinder of 
finite spatial size (circumference) L. Bethe equations become infinite in number in this 
limit. They must be substituted by some sort of "density of roots" method. However the 
traditional density approach, based on linear integral equations, is valid only in the large 
L limit and does not take into account finite size effects. The NLIE approach, instead, is 
able to fully control the physics on the continuum at any size L. 

The lattice NLIE (j3J) has a sensible scaling limit only if we fine tune the inhomogeneity 
parameter 9 to go to infinity as In N in a very precise way 



6 - In 



2N 



thus introducing a mass scale M.. The only variation in the NLIE is in the Dj, term that 
becomes simply 

D b (9) = MLsmhO 

The asymptotic particle interpretation of this term will be clear in the IR analysis of 
next section. We are in presence of a genuine renormalization procedure, leading to a 
renormalized quantum field theory on the continuum. In the following it is convenient to 
introduce the dimensionless scale parameter £ = M.L. The UV physics of the model is 
reproduced for t — > (negligible mass scale), while the IR one corresponds to the "large 
volume" £ — » 00 regime. 

The counting equations on the continuum are slightly simplified with respect to the 
lattice. They read 



N 



H 



2N S = M C + 2(S + M W + 

= S + M W + M SC + M [ c ] -M R - 5 



(22) 



(23) 
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where S is now the topological charge of the continuum theory, connected to the lattice 
spin S z is: 

2S = 2S Z - N+- AL. (24) 

The terms with integer parts on the rhs. of (J2*UJ) present on the lattice and coming from the 
value of Iml?(±oo) become zero in the continuum, where B(±oo) = 1. The same happens 
in the additional parity constraint for Ni that now simply becomes the requirement that 
it must be even: Ni G 2Z. 

The NLIE can be solved numerically with an iterative procedure for a given set of 
quantum numbers Once the position of the sources and the functions ln6, \ny have 
been calculated, the total energy and momentum of a state can be calculated through the 
formulae 



N H n h M S 

E = M ( cosh(hj) - 2 ^ cosh(/if ) - ^ {cosh( Sj ) + cosh(s~ )} (25) 
, j=i i=i i=i 

+oo +oo 



M c . . . . 

Y cosh( Ci ) + — J d6smh(e + ie) In B(6) - J d6 sinh(6 - ie) \nB(0) 



-oo — oo 



' N H N h m s 

P = M [Yl s ™h(hj) -2^2 sinh(/i?) - ^ {sinh( Si ) + sinh(s})} (26) 
, i=i i=i j=i 

+oo +oo 



- V sinh(c,) + — / d0cosh(0 + ie) In B{B) - — / ci0cosh(0 - ie) In £(0) 

~^ 27T ,/ 27T J 



Notice that wide and self-conjugate root contributions drop from this expressions because 
sinh// 6 = cosh// = in the repulsive regime p > 1 we are considering here. 



5 Large scale analysis and Particle Theory 

It is interesting to verify which scattering data can be reproduced by the large £ limit of 
this construction. In the £ — > oo limit the convolution integrals involving B and B in both 
equations are exponentially depressed, as can be easily seen by replacing the leading large £ 
contribution that goes as £sinh0 into the convolution itself. However, this same reasoning 
cannot be applied to the other convolution integrals involving Y, as the function y has no 
exponentially leading term. As a result, the I — > oo asymptotics of the NLIE looks like 

-% In 6(0) = C b + £ sinh + gi {9) + g h {B) - i{K [ ~^ +e] * In Y){9) + 0(e~ l ) (27) 
\ny(6) = C y + g y (6) + 0(e- e ) (28) 
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The dominating term of the first equation for large £ is £sinh0, thus showing monotonicity 
of the real part of the function —i In b, which in turn implies that for large I there are no 
special objects: N s = 0. 

The second equation is now purely algebraic and can be solved for y and the result can 
be plugged into the first equation. Defining 

K{9) = (K^ * ln(l + e c » +a *))(9) 

we have, putting u = and letting 5 b take both its values or 1 

-i In b{9) = £ sinh 9 + g b (9) - iK(6) + 5 b 7r (29) 

The energy and momentum expressions also simplify in the same way 

N H M c 

E = MJ2 cosh(^) - M cosh(c,-) + 0(e~ e ) (30) 

3=1 3=1 



N H M c 

P = MJ2 sinh(^) - M ^ sinh( Cj ) + O(e^) (31) 

3=1 3=1 

Let us first examine the case of Nh holes without any complex pair (other than, of 
course, the quasi 2-strings of the vacuum sea, that do not contribute here, being confined 
in the dropping convolution integrals). The energy and momentum are given by 

N H N h 

E = M^ cosh(hj) , P = M Y smh(hj) 

3=1 3=1 

We can interpret this result by saying that holes represent a system of asymptotic particles 
of mass M. and rapidity hj. 

The quantization of particle rapidities for large volume in a factorized scattering theory 
is given by the formula 

e iMiainh9k A{9 k \{9i}) = 1. (32) 

where A(0|{0j}) is an eigenvalue of the n-particle color transfer matrix defined by the 
products of two-body S-matrices 



T(9\{9 i })=l[S(9-9 i 



i=i 



A comparison of (|32|) with the exponentiated form of (|29|) and recalling the quantization 
conditions {THJ...IIEJ) shows that in general 

A(9 k \{9, l }) = e i9b ^ ) e lc{ek \-l) Sb 
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In the case of two particles, A are simply the eigenvalues of the two-body S-matrix. 
Moreover, if the S-matrix has a tensor product form, also its eigenvalues have. The two 
particle states are realized by choosing N H = 2. The two hole positions are denoted hi,h 2 . 
We recall the well known fact that 

= S (9) 

where S (9) is the S++ = SZZ element of the SG two-body S-matrix. The counting 
equation (|22j) gives a few possible choices 

• S = 1 and M c = M w = M sc = 

• S = and M c = 2, M w = M sc = 

• S = and M sc = 1, M c = M w = 

Notice that Mw and Mc can only be even, as the complex roots other than self-conjugate 
always come in pairs. From the other counting equation (|2~3l one can see that for all the 
2-hole states N t = 0. 

A simple inspection case by case of the first factor e tgb ^ k ^ of A shows that it reproduces 
exactly the SG S-matrix and the quantization of close pair imaginary parts studied in jE]. 
We conclude that 

\ = \ SG e iK {-lt 

The Ahn S-matrix will be reproduced if e lK [— Vf b gives the four eigenvalues of the RSOS(2) 
S-matrix. As 

y (9) = (-l) 5 y tanh I J tanh \—-* 

for the calculation of K we have to distinguish the two cases 5 y = and 5 y = 1. 
In the 5 y = case one gets 

K{B) = iTi {Q 1 (9 - h 12 ) - Q 2 (9 - h) - Q 2 {9 - h 2 )} , (33) 

where h 12 = hl ^ h2 and 

1 i 
QJ0) = -— arctansinh(^) - — lncosh(^), (34) 

where X2(9) is equal to x(9) with p = 2. Important values from the point of view of the 
quantization equations are 

K{h l ) = n{h 1 -h 2 ), (36) 
K(h 2 ) = Tl{h 2 - h), (37) 
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where 

n{B) = ~X2(6) - '-MO), (38) 

sinh^ , 
= (39) 

The following important identity yields 

Xa{6) = \x2{6) + \m0)- (40) 

where x±{Q) denotes x{@) calculated at p — 4. This gives two of the four eigenvalues of the 
RSOS(2) S-matrix, namely 

. sinh(^) 

sinh(H^) 1 ' 

Along similar lines, the case 5 y = 1 gives the other two eigenvalues 




thus exhausting the RSOS part of the two-body S-matrix. The complete structure of the 
solitonic S-matrix of Ahn [TH] is reproduced. 

We observe that the eigenvalues (I41|42l) of the RSOS part of the color transfer matrix 
do not contain the 2 e ^ 2m crossing factors of Ahn's S-matrix. This is because the crossing 
factors correspond to a gauge transformation of the RSOS transfer matrix, under which 
the eigenvalues obtainable by finite size effects are invariant [33] . 

The main surprise of this discussion is the fact that to take into account all the 4 
RSOS eigenvalues we had to take into account both values 5 y = 0, 1. Now, from the lattice 
construction, for the states here examined, only S y = is allowed. The case where in the 
second NLIE 5 y is taken to be 1 cannot be deduced from lattice construction. To say that 
there exists a NLIE where 5 y = 1 is a conjecture, similar in a certain sense to the one 
leading to the NLIE for odd topological charge sectors in SG jHj- Allowing this extension 
of the NLIE does not only reproduce the two missing eigenvalues of RSOS part of the 
S-matrix, but, as we shall see in the following, also the entire R sector of the UV modular 
invariant partition function. If there is a viable lattice construction for this choice of 5 y or 
not remains to be clarified. 



6 UV limit and kink approximation 

Now we examine the ultraviolet limit of the states described by the NLIE. We follow the 
standard approach described in detail in [H [20111]. The position of the sources for I —>■ 
can remain finite (central objects), or they can move towards the two infinities as ±ln (|) 
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(left/right movers). We introduce the finite parts 0- ' of their positions 9j by subtracting 
the divergent contribution: 

W-{^±lnQ),^°}. 

We denote the number of right/left moving and central objects by N^'°, Np°, M c '°, . . . 
etc. For later convenience we introduce the right/left moving and central spin given by 

S±>° = l - {Np° - 2Np° - M±'° - 2M±>° - 2Mff) . (43) 

In the UV limit the NLIE splits into three separate equations corresponding to the three 
asymptotic regions. This is why for all the auxiliary functions of the NLIE (jlj) we define 
the so called kink functions as 

F±(6) =JirnF \9±\n^j . F e {log 6, logy, logo, . . . }. (44) 

In the UV limit these kink functions satisfy the so called kink equations and the energy 
and momentum can be expressed by them. Performing the above kink limit on our NLIE 



the kink equations take the form 

log&±(0) = iCll±ie ± V +ie ) + ig^{9 + ie) + ig®(9 + ie) 

+ (G*lnB ± )(6)-(G +2e *lnB ± )(0) + (K-% +e *\nY±)(9), (45) 

logy± {9 - i|) = C y± + gfl(9) + * \nB ± )(9) - {K^ * \nB ± )(9), (46) 

logj/±(0) = lim logy± [9 - £ + irA . (47) 

lna±(0) = C- a + iC± + i~g < il{9) + (G7/ * lnB ± )(0) - (G+/ * In B ± ){9), (48) 



where in (|4*5|l i G {1,2} depending on the value of the imaginary part of 0. If < Im0 < 
7r/2, then i = 1, if tt/2 < Im0 < ir, then i = 2. The bulky expressions of the constants 
and the concrete form of the source functions are collected in appendix A. The finite part 
of the right and left moving objects can be obtained from the kink functions by imposing 
quantization conditions very similarly to ()13| ... fT8|) . 

After some usual manipulations [H [201 HI it turns out that in the UV limit the energy 
and momentum can be expressed by a sum of dilogarithm functions with the — > ±00 
limiting values of the kink functions in their argument. One of these two limiting values is 
trivial, namely: 

&±(±oo) = &±(±oo) = 0, (49) 
j/±(±oo) = (-I) 5 ". (50) 
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The other limiting values satisfy a nontrivial coupled nonlinear algebraic equations, the 
so-called plateau equations [T]: 

log b ± (Too) = iC b ± 2i Xoo (S - 2S ,± ) ± 2m Z± ± in (m sc + ^--Nf 

+ ^ {ln5 ± (Too)-lnS ± ( T oo)} + Jlny ± (Too), (51) 

2/±(Too) = (-1) 2S± +^ (5 ± ( T oo)£ ± ( T oo))*, (52) 

where 

M ± 

$ = &--f- Mf c , (53) 

irrelevant from the point of view of the exponent of equation (|51jl. Now we have to solve 
eqs. ()51|52jl analytically. Motivated by the TBA formulation of the problem, we take the 
following Ansatz for the solutions of the plateau equations (|51l52j) : 

b±(Too) = e ±3ip± 2 cos(p±), &±(=Foo) = e T3ip± 2 cos(p±), (54) 



5±(Too) 



D ±2i P± sin (3p±) 
sin(p±) ' 

sin (3p±) 



B±(Toc) 



sin (3pj 
sin (p± 



y ± ( T oo) = 4cos(p±) 2 > 0. 



(55) 
(56) 



sin(p±) ' 

Since we need only &±(=Foo) and not the (extended) log of it, we can restrict the allowed 
values of p± in the [0, 2tt] interval. The solutions of the plateau equations formally take 
the form 

p ± = vr (kf ± 5 b ± - + A p± ) - (2k± ± 25 b ± 2- + 3A 



where 
A 



S - 2S + - Nt, 



N: 



± 



3^ - 


P± 


n i = 


2^ 




P± 


- 7T - 


- 71 - 




- 71 - 




- 71 - 



and further constraints must be satisfied by the parity of the integers ki^ and N^: 



(57) 
(58) 

(59) 
(60) 

Using the dilogarithm sum rule of appendix B. and putting everything together the con- 
formal weights take the form 



fcj = 2ZJ + M sc + ^-Nt 



n 



± 



Nf mod 2 = 2S ± + 8 y mod 2 



16 y 2 V 2 



AT. 



(61) 
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where 



n ± = -[5 b ±kf±^A p± ), (62) 
2^- -p iy p 



N± = Np — =p S ± n± - l(S - S ± )S ± J N± = N± mod 2. (63) 



J± = ±4± =f 2J fc ± =f 2J,± =f 2/ c±T =f 2I tl) ±t =F l w ±j ± 4(D± 

?± i Jl/fi i I\,f±\ / Jf _ 1— -r- J\r± _L /l/r -i- Jl/f± 



T (S* + M± + Mf c ) I8 b ± -± =p iVf ± M sc =F M± j (64) 
- (M++M+) (S-S^+i^. 
where if + and iT_ are given by 

K + = k++ Mg )+ (JVl - iV+) - ^- {(5, + 2(5 + M w + M sc ) + - M { J )+ } (65) 



K- = kw + ~f (*„ ~ Mg ] + M®" ) (66) 

where depend only on the actual configuration of wide roots, and they are integers, if 
Mf c are even, and they are half-integers if Mf c are odd. 

Analyzing the formulae (|61H66|) of the UV conformal weights, one can see that S can be 
identified with the winding number m, the parameter 5 y G {0, 1} of the NLIE distinguishes 
the NS and R sectors of the theory, and depending on the state under consideration the 
sum N± + J± can be either integer or half-integer. Furthermore it can be easily proven 
that there is a relation between n± of |62l) and the winding number S, namely in the NS 
sector 

n±EZ if S E 2Z, (67) 
n ± GZ + ^ if SE2Z + 1, (68) 

while in the Ramond sector 

n±eZ + ^ if S G 2Z, (69) 



n± G Z if SG2Z+1. (70) 

in full accordance with what described on the UV operator content in section [2j Moreover 
it can be proven that in the S y — 1 case (R sector) the sum N± + J± is always integer 
as it must be in the R sector of a c = | CFT. Thus, due to the previous remark and the 
relations (|67H7()L the conformal weights (|61H66jl can be interpreted within the framework 
of c = | CFT as conformal weights appearing in the modular invariant partition function 
of Di Francesco et al. [TT| . 
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7 Comparison with PCFT results 



In the previous section we saw that the NLIE (jlj) describes a modular invariant c = | 
CFT in the UV, but we would like to have more evidence that our coupled NLIEs describe 
the SSG theory. The simplest check is the numerical comparison of NLIE and perturbed 
conformal field theory for certain states. This comparison is really important in the R 
sector where the form of the NLIE is a pure conjecture and can not be obtained from the 
integrable light-cone lattice regularization of the model. 

As we already discussed in sectional the SSG model can be formulated as a perturbed 
CFT with the action on the cylinder 

A = A 3/2 + g J d 2 z $ pert (z,z), (71) 

where A3/2 is the action of the c = | CFT, g = is the coupling constant with mass 
dimension y = 2(1 — l/R 2 ) with R being the compactification radius of the conformal 
normalized boson field and the & pert (z,z) perturbing operator of eq.j2} is normalized by 
(Q P ert(z, z)$ pert (0, 0)) = \z\~ 2( - 2 ~ yS) . In the framework of conformal perturbation theory the 
energy levels can be represented as a sum of integrals over the n-point functions of the 
unperturbed CFT [22]. In the case of SSG theory there is no bulk energy contribution on 
account of supersymmetry, and the perturbation series contains only even powers of the 
coupling g. Using the exact relation between g and the massgap of the model [2 



9* P 2 \ ,,i f> 2 (3 2 V ^ , , T(x) . , 



the perturbation series can be expressed as power series of £ 2y , where t = M.L is the 
dimensionless volume. The first correction to the energy of a state \a) having conformal 
dimensions A a , A a can be written as 

— E a (£) = -(- - 12(A a + A a )) - C l 2 a) fy + 0(£), (73) 

7T I 

where the leading order coefficient reads as 

Cf = a— (74) 

7T 

where 

contains the massgap formula and is a simple integral expression 

f d 2 z 

I la) = J T^(a\%ert(l,l)%ert(z,z)\a) conn . (76) 
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In (THJ the first nontrivial PCFT corrections were calculated for various states in the NS 
sector and they were used to test the TCSA near UV, and to check the NLIEs proposed 
to describe the three ground states of the model. Now, having NLIE for excited states as 
well, it is worthwhile to check it against PCFT predictions. Unfortunately the comparison 
is quite restricted on both sides. Actually, the applicability of PCFT, mostly in the R 
sector, is strongly restricted by IR divergencies. Thus, up to two-kink states the number 
of IR safe states is very small. On the other hand the NLIE can not be solved numerically 
with present techniques if special objects appear [T]. Thus, the PCFT-NLIE comparison 
can be done only for a few states in certain restricted regions of the coupling constants. 
In the NS sector, in the zero momentum two-kink states space, the energy of the state 

|l-> = ^ (y§°\0,0)-V±°$(0,0) \ |0) with conformal weights A = A = ^ can be 
calculated numerically from NLIE in the region 1 < p < 2. The PCFT correction of the 
state was calculated in [18j and is given by 



^ = 7 ( ; 



2 2R 2 



1 \ (\ 3 \ 1/1 1 \ / 1 



R 2 1 \ 2 2R 2 2 7 \2 2R 2 ) 1 \R 2 



(77) 



In the language of Bethe Ansatz this state corresponds to a S = state with 2-holes, one 
self-conjugated root, and 5b = 5 y = 0. We checked numerically our NLIE against PCFT 
at a lot of values of the compactification radius R, and in every case we experienced very 
good agreement. To illustrate the agreement between NLIE and PCFT, the numerical 
comparison of the quantity e(l) = ^E a (l) + (| — 12(A a + A a )) coming from NLIE and 
PCFT at the specific R 2 = | point can be found in table 1. One can see that the two sets 
of data approach one another as the volume tends to zero. 

In the space of zero-momentum two-kink states of the Ramond sector, because of IR 
divergencies, the only state for which PCFT can be applied is a charged 2-kink state. At 
UV, this state corresponds to \S = 1,R) = -Ro,i(0, 0)|0) with conformal weights A = A = 
^ + | (f ) ) an d the first nontrivial correction to this state is given by 

^* = _- 7 (_i + ^) 7 g + ^) 7 (__y. (78) 

In the language of Bethe Ansatz, this state corresponds to a pure two-hole state with 
5b = 0,5 y = 1. The numerical comparison of e(£) coming from NLIE and PCFT at the 
specific value R 2 = | can be found in table 2. The two sets of data converge as the volume 
tends to zero. 

The very good agreement of NLIE and PCFT data give an additional strong support 
that the NLIE (J3J indeed describes the finite size effects of the SSG model. 



8 Conclusions and Perspectives 

In the present paper we have presented the full excited states NLIE for the N=l super- 
sine-Gordon model, deducing it from a lattice NLIE introduced by Suzuki for the spin 1 
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£ 


e{£) (NLIE) 


e(£) (PCFT) 




0.1 


0.213251889(1) 


0.2135198827 


5- 10~ 3 


0.05 


0.096639933(5) 


0.096695005 


lO" 3 


0.01 


0.015365337(5) 


0.01536672004 


2 ■ 10~ 5 


0.005 


0.0069587125(5) 


0.006959000 


5 • 10" 6 


0.001 


0.001105913(5) 


0.00110592082 


io- v 


0.0005 


0.00500832(5) 


0.005008293 


3- 10~ 8 


0.0001 


0.000079592(5) 


0.00007959153 


7-10- 10 



Table 1: Numerical comparison of PCFT with NLIE for the state |1-) at R 2 



i 


e{£) (NLIE) 


e{£) (PCFT) 


ty 


0.1 


0.4411699(1) 


0.42645819 


0.025 


0.05 


0.24978747(1) 


0.24493591 


8 • 10~ 3 


0.01 


0.06795829(1) 


0.067589069 


6 • 10" 4 


0.005 


0.03894153(2) 


0.03881972 


2 • 10" 4 


0.001 


0.01072140(3) 


0.010712145 


10~ 5 


0.0005 


0.00615559(2) 


0.00615251 


5 • 10~ 6 


0.0001 


0.00016979(1) 


0.000169776 


4- 10- 7 



Table 2: Numerical comparison of PCFT with NLIE for the state \S=1,R) at R 2 
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integrable XXZ chain [Hj. We have checked the consistency of this procedure by comput- 
ing some of the infrared and ultraviolet properties known for such quantum field theory. In 
particular, the S-matrix proposed by C. Ahn [TH] has been reproduced, as well as various 
properties of the c = § CFT that is expected to describe the ultraviolet behavior of the 
model. Much more work has to be done, however. First of all, contrary to the sine-Gordon 
case, we lack here of a deduction of the equation of motion from the lattice construction. 
This keeps the whole construction quite conjectural (although on a very solid ground). A 
neat deduction of equations of motion from the light-cone lattice construction or alterna- 
tively from a bosonization of the alternating inhomogeneous XXZ Hamiltonian could be a 
very important step forward. 

We have investigated in detail only the solitonic sector of the model, by focalizing on 
the repulsive regime where no breathers appear. The structure of the supersymmetric 
breathers should be investigated for at least two reasons: it completes the analysis of the 
asymptotic particle spectrum of the theory and gives access to quantum group restrictions 
describing superconformal minimal models perturbed by $13. 

The N=l SSG model is just the first in a whole series of Fractional SSG models [2H] that 
can be seen as perturbations of the conformal models introduced in [TJj. The full series 
of coupled NLIE's generalizing those presented here has been proposed, for the vacuum 
state, by C. Dunning [22] . These models are evidently mapped, by analogy with SG 
and SSG, into light-cone higher vertex models like the 44-vertex, etc... or equivalently in 
alternating inhomogeneous XXZ chains of higher and higher spin. It would be interesting 
to check Dunning's conjecture and extend NLIEs to the structure of excited states. Two 
applications of this could be of importance: the quantum group restriction leading to 
S77(2)-coset theories (SH] perturbed by (p^l and the large spin limit where this series should 
make contact with the N=2 SSG model, whose finite size effects are of principal importance 
in the formulation of superstrings propagating in pp-wave background pTT| . 

An even more unifying point of view has been proposed in [2H], where two massive 
coupled NLIEs encompass the whole set of SU (2)-cosets, of FSSG models and their parent 
theory, the so-called SS-model [32J. 

Finally an investigation of the SSG model with integrable boundaries, along the lines 
of what done for SG in (2H1 would be of interest. 

We hope to return to these issues in future. 
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Appendix A 

This appendix is devoted to list the constants and source functions appearing in the kink 
equations (I45II48I) . The list of values of the constants is as follows 



eg = C b + 2 X oo (S -S ± )±n (M w - M± + M sc - M±) + 2vrZ± ± |(iVi - Nf), (79) 

C$ = C$-2ir{N x -Ntl Cf^cff, (80) 
(%+ = z 7T + 2(5 + M w + M sc ) + Mg ] - (fi + + M+ + M+ + M ( J )+ ) } , (81) 
<V = tTT {5, - M^ 2) + (fir + M w + M" + M®")} . (82) 

C7 g = i|2wU + ~J + ^(iV + -AL)}. (83) 

C~ a± = ±4 Xoo (fi - fi±) =f 27r(fi - fi±) + 2^, (84) 

The numbers and l w are integers which appear only in the presence of wide roots, and 
their actual value depends on the relative position of the wide roots. The form of the 
source functions of (|45H48j) are given by 

&W = Y,Xx(0-h?* fl ), (85) 

E Xk(6 - hf) - 2 E - ^ -EW»- #°) + Xk(^ - 
i=i i=i i=i 



j=i j=i j=i 

M ±fi M ±fi M±>° 

-_c l u w III. 

E X(0 - of) - E Xn(0 - - E Xn(9 - ^ c )± '°), (87) 
j=i j=i j=i 
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A' 



±.o 



40) 



±,0 

E x»( 



N 



S±,0 



Mi 



±,0\ 



±,0\ 
C i ) 



2 E xn( 



S±,0\ 



E 



±,0n 



-±,0 



i=i j=i 



w 



(i)±,o 

sc /J // 



)) 

(88) 



Appendix B 



In section in the calculation of the UV conformal weights, the following dilogarithmic 
sum must be calculated: 



S (p) =2\L + [e iip 2 cos(p)] + L + [e~ 3tp 2 cos(p)] + L, 



sin (3p) 
sin (p) 



where 



£+0*0 = \ J dy | 



ln(l + lny 



2/ i + y . 

The function S'o(p) has the following simple properties: 

S (p) = S (p + 7r), S (p) = S (-p). 
After some algebra one can prove that 

iV p 7r In (4 cos 2 p) , 



2tt 2 

S (p) = — + N„l2n 



p — IX 



where 



N p = N p mod 2 7V P 



[3^1 




"P" 


L 7T- 




-7T- 



(89) 
(90) 

(91) 

(92) 
(93) 



and we made the choice of ln(— 1) = in. 

Appendix C 

We solve the Bethe ansatz equation on the lattice directly for N = 36, 7 = ^ and for 
various particular £ = 2Ne~ e in the ground state. This corresponds to the attractive 
regime and 5 y = 0. The Bethe ansatz roots and ^sm(\nb) are plotted below. Here we use 
the variable x = jO/n. 

One sees that the central regime in the graph of In b becomes flatter for smaller £. The 
slope however seems to be always positive. 

With increase of N, the Newton method for finding the Bethe ansatz roots becomes 
quite unstable. We thus adopt a kind of WKB method to locate them. We plot the result 
for N = 128 in Fig. El This again support the positivity of the slope qualitatively. 
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Figure 1: The Bethe ansatz roots (left) and the corresponding 3m (In 6) for N = 36,7 = 
tt/2.5,£= 1 




25 



Figure 3: The Bethe ansatz roots (left) and the corresponding 9R(ln&) for N = 128,7 = 
7r/2.5,^ = 0.1 




Figure 4: The Bethe ansatz roots (left) and the corresponding 3m (In 6) for N = 38,7 — 
tt/2.5,^ = 0.1 
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Figure 5: The Bethe ansatz roots (left) and the corresponding (In b) for N = 38,7 = 
7r/3.5,^ = 0.1 




Next we solve the BAE for N = 38. This corresponds to the choice 5f, = 1. The plots 
in FiglUshow that the slope near the origin is again positive (modulo jumps) for 7 = 7r/2.5 
for smaller values of I. 

The plots in the repulsive regime (7 = 7r/3.5) are given in Fig. El The right figure shows 
the negative slope at the origin. (Imagine the "extended Brilliouin zone" and line defined 
by Smi = 7T.) This indicates that the upper component 6j of the quasi 2-string centered 
at the origin has j|lmln&(0)| < 0, i.e. it is a special object. While in the analog 
attractive case N s = and for the vacuum where there are no holes N e ff = N H — 2N S = 0, 
here we have now N s = 1. However, the jumps at x ~ ±0.7 do not correspond to any pair 
of roots. Therefore they identify two pairs of holes (they are not special, as the derivative 
of In b is positive for them). In total we have now Nh = 2, but N e ff continues to be 0. 
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